I. COLLECTIVE EXCITATION SPECTRUM OF DENSITY MODULATED ONE 
DIMENSIONAL ELECTRON GAS (DM1DEG) 



In the last decade or so remarkable progress has been made in epitaxial crystal growth 
techniques which have made possible the fabrication of novel semiconductor heterostruc- 
tures. These modern microstructuring techniques can be used to laterally confine quasi-two- 
dimensional electron gas (2DEG) in e.g.,GaAs/AlGaAs heterostructure on a submicrometer 
scale to quasi-one dimensional structures (quantum wires) or quasi-zero-dimensional quan- 
tum dots. The system that we are considering in the present work can be realized by 
various methods, e.g., the application of a laterally microstructured gate electrode and the 
holographic illumination technique which allow for a tunable periodic density modulation of 
the quasi-one-dimensional electron gas. 

There have been a number of approaches to study collective excitations s pec trum of a 
quasi-two and one dimensional electron gas (1DEG) systems theoretically 
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and experimentally [12|, Il3l ? ] . We extend the theoretical work by calculating the inter- and 
intra Landau band magnetoplasmon spectrum of a density modulated quasi-one dimensional 
electron gas (DM1DEG) in the presence of a perpendicular magnetic field using the Self- 
Consistent-Field (SCF) approach and focus on the oscillatory behavior of the intra-Landau 
band magnetoplasmons. In this context, the term "ID" means that we start with the original 
2DEG in x—y plane. We apply a confining potential in the x-direction leaving the y-direction 

free. A magnetic field is applied along the ^-direction perpendicular to the x — y plane of 

h 2 k 2 

the original 2DEG. The collective particle energy spectrum, Eij(k y ) = 2— if + E l x consists 
of energetically separated ID subbands formed due to lateral confining potential along the 
x-direction. The electron wave vector k y characterizes the free motion in the ^-direction. 

The effect of density modulation is to broaden the Landau levels into minibands whose 
width oscillates as a function of the magnetic field strength. The electronic states are thus 
substantially altered, resulting in modulated density of states, as shown by magnetocapaci- 
tance measurements [12I, [3] of the quasi-two-dimensional systems. Behavior akin to this is 
expected for the quasi-one dimensional system under consideration. The density of states 
affects many response and transport phenomena as well as thermodynamic properties. Of 
these, one of the most important properties is the collective excitation spectrum and we 
evaluate the dynamic, nonlocal dielectric response function to study it. As we show in our 
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work for a quasi-one-dimensional system, magnetoplasma spectrum significantly exhibits 
modulation of the electronic density-of-states as oscillating magnetoplasma frequencies. As 
we discuss in detail, this result is obtained in the regime of weak modulation and long wave- 



length thus we don't find this result discussed in [26( . Another condition for the observation 
of these oscillations is that the coupling between intra-Landau band mode and inter-Landau 
band mode must be small. Mixing of these modes can be minimized by controlling the de- 
gree of density modulation and by applying an appropriate magnetic field and confinement 
potential. 



II. FORMULATION 

Our system is a density modulated quasi-One Dimensional Electron Gas (DM1DEG) 
built on a two dimensional electron gas (2DEG) by inducing another effective confining 
potential \(m*uj 2 x 2 ) along the x-direction which is assumed to be parabolic. The magnetic 
field is perpendicular to the x — y plane in which electrons with unmodulated areal density 
rip, effective mass m* and charge — e are confined. We employ the Landau gauge and write 
the vector potential as A = (0,Bx,0). The two-dimensional Schrodinger equation with 
parabolic confining potential in the Landau gauge is (h = c = 1 here); 

ff »=it-^ + H| + >f] + i(^oV), (i) 

Since the Hamiltonian does not depend on the y coordinate, the unperturbed wavefunctions 
are plane waves in the ^/-direction. This allows us to write for the wavefunctions, 

4> nky (x) = -^=e lk « y u nky{x) , (2) 
V L y 

with L y being a normalization length in y-direction and x a 2D position vector on the x-y 
plane. Substitution of the above form of the wavefunction in to equation (1), yields 

with (fljj = u 2 + u 2 ,l 2 (B) = l 2 ffi), where cu c = ^ is the cyclotron frequency, x = 



—l 2 {B)k y = — ; is the coordinate of cyclotron orbit center and 1(B) is the mag- 
netic length, m*(B) = rn*^ is the normalized effective mass. In the x-direction, the 
Hamiltonian has the form of a harmonic oscillator Hamiltonian. Hence we can write 
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the unmodulated eigenstates in the form (j) nky (x) = -j=e tkyV u nky{x . XQ) , with u nky{x . Xo) = 
(^2 n n\l)^ L exp(— 272(0; — XQ) 2 ~)H n {^—f^) , where u nky ( x - xo ) is a normalized harmonic oscilla- 
tor wavefunction centered at Xq and H n (x) are Hermite polynomials [15] with n the Landau 
level quantum number. The energy of nth Landau level (unperturbed Hamiltonian Hq) is 

s m = ( B + i/2 )no + 5 |_. (4, 

In the presence of modulation, the Hamiltonian is augmented by the term H' = Vq cos(^x), 
where Vq is the amplitude of the modulation and is about an order of magnitude smaller than 
Fermi energy {Vq/ef)- Due to smallness of Vq we employ first order (in H') perturbation 
theory in the evaluation of the energy eigenvalues using the unperturbed wavef unctions. The 
correction to the unperturbed eigenenergies with a new variable y = is given by: 



^—cos(— x Q ) [ dy exp(-y 2 )[H n (y)] 2 cos(— ly). 
n n\ a J a 



V^2 



The integral in above equation is given by Gradshteyn-Ryzhik [16] (page 841 # 7.388.5); 
and the result is; 

= V n cos(^-x ), (5) 

where V n = V exp(— X/2)L n (X), X = (^ L ) 2 2r ^ e ^ ,and L n (x) is a Laguerre polynomial. By 
combining equations (4) and (5) we write, for the energy eigenvalues to first order in H', 

k 2 2tt 

e{n, x ) = (n + 1/2)Q + - y + V n cos(—x ). (6) 

2m* (B) a 

The above equation shows that the formerly sharp Landau levels, equation (4), are now 
broadened into minibands by the modulation potential. Furthermore, the Landau bandwidth 



V n |) oscillate as a function of n, since L n (X) is an oscillatory function of its index 15] 



III. DENSITY-DENSITY CORRELATION FUNCTION OF A DM1DEG IN A 
MAGNETIC FIELD 

The dynamic and static response properties of an electron system are all embodied in the 
structure of the density-density correlation function. We employ the Ehrenreich- Cohen Self- 



Consistent Field (SCF) approach J9] to calculate the density-density correlation function. 
The SCF treatment presented here is by its nature a high density approximation which has 
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been successful in the study of collective excitations in lower- dimensional systems such as 
Planar se— ctor superb Q and q _ whe st_ 0, flfi, both with 
and without an applied magnetic field. Such success has been convincingly attested by the 
excellent agreement of SCF predictions of plasmon spectra with experiments. 

Following the SCF approach, the density response of electrons due to a perturbing po- 
tential is given by 

X I- 4\ f( £ a') — f(£a) ^ \\r(- \\ 

on(x , z :t) = > < a \ V (x, z:u>) \ a > 

^ So; - e a + u + irj 

aa' 

x < a' I 5(x — x Q )S(z — zq) I a >, (7) 

where V(x, z\ uj) is the self-consistent potential and a stands for the quantum numbers n 
and k y . Fourier transforming on the x — y plane we obtain the induced particle density 

/(£<*') - /(e«) 



5n{q, z ; u) = -S(z)V(q, z = 0; uj)^ Z 



A ,£ a > - e a + uj + ir) 

aa' 



x |< a' | e ^ | a >| 2 , 

where A denotes the area of the system. We can perform the fc^-sum in the above equation 
to obtain 

6n(q, z; uj) = —5(z)V(q, z = 0;uj) 

V r ( , liE^Lih^y)} ~/( £ ( w > fc t/)) 

^ nn ' [ 2m*n 2 ) e(n' J k y -q y )-e{n J k y )+uj + z V {) 

n,n' ,k y " " " 

Writing the induced particle density as Sn(q, z; uj) = Sn(q, u)S(z), allows us to rewrite equa- 
tion (8) as Sn(q, uj) = V(q,uj)Yl (q,uj): where V(q,uj) = V(q,z = 0;uj) and Ho{q, oj) is the 
density-density correlation function of the non-interacting electron system, given by 

n (n < ,\ - -X^X^r < ^ c ^ f( £ ( n '> k v-%)) ~ f^(n,k y )) 

n oW , uj) - A 2^^ nn '^2m*nl } e{n', k y - q y ) - e{n, U+u + i-q' 1 ' 

where 

C nn ,{x) = ^e- x x n ^[L n n \~ n ^x)f 
n\\ 

with rt\= max(n, n'), n,2= min(n, n'), and L' n (x) an associated Laguerre polynomial. The 
induced potential V md is related to the density response Sn by Poisson's equation 

V 2 V ind (x, z; t) = — 5n(x, z; t), (10) 



where k is the background dielectric constant. The above equation can be solved to yield 

Ittp 2 

V ind (q,u) = —5n(q,u;). (11) 

Recalling that the self-consistent potential, V(q, u) = V ext (q,u) + V md (q, u), is the sum 
of the external and induced potentials, multiplying both sides by H(q,u) and solving for 
6n(q,u) yields 

Sn(q,u)=U(q,u)V ext (q,u), (12) 

where 

W ' ; l-v c {q)U {q,u) y J 

is the density-density correlation function of the interacting system with v c (q) = the 
2-D Coulomb potential. Making use of the transformation k y — > — k y with the fact that 
e(n, k y ) is an even function of k y , and at the same time interchanging n <-> n' we write for 
the non-interacting density-density correlation function equation (8) 

a 

U (q,uj) = -y]C nn /( ° ) / dx [f(e(n,x + x' Q ) - f(e(n',x ))} 

o 

x [e(n, x + x' Q ) - e(n', x ) + uj + it]]' 1 . (14) 

In writing the above equation we converted the ^-sum into an integral over xq. f(E) is the 
Fermi -Dirac distribution function, xq = — ky ^ and x' n = — 9v ^ . 

The above equations (13,14) will be the starting point of our examination of the inter- and 
intra-Landau band plasmons. The form of the expressions for the real and imaginary part of 
the density-density correlation function makes the even and odd (in frequency ui) properties 
of these functions very apparent. These functions are the essential ingredients for theoretical 
considerations of such diverse problems as high frequency and steady state transport, static 
and dynamic screening and correlation phenomena. 

The plasma modes are readily furnished by the singularities of the function IT(g, a;), from 
the roots of the longitudinal plasmon dispersion relation obtained from equation (13) as 

l-v c (q)ReU (q,oj) =0 (15) 

along with the condition /mn (g,w) = to ensure long-lived excitations. The roots of 
equation (15) give the plasma modes of the system. 

2ne 2 m*^}^^^^ ' ^ 2UJC 

n,n' 



i^£ c -(^>('H +'<-«)>■ < 16 > 



with 

a 

7Y,,^ = P / /*rv, 

+cu - e(n, x ) + e(n', x + £q) 
o 

where P is the principal value. 

The plasma modes originate from two kinds of electronic transitions, those involving 
different Landau bands (inter-Landau band plasmons) and those within a single Landau- 
band (intra-Landau band plasmons). Inter-Landau band plasmons involve the local ID 



magnetoplasma mode and the Bernstein-like plasma resonances [20], all of which involve 
excitation frequencies greater than the Landau-band separation (~ Oo) ■ On the other hand, 
intra-Landau band plasmons resonate at frequencies comparable to the bandwidths, and the 
existence of this new class of modes is due to finite width of the Landau levels. These mag- 
netoplasmons in quasi-one-dimensional system have been analysed in detail elsewhere 26]. 
We will concentrate on the oscillatory behavior of these magnetoplasmons. The occurrence 
of such intra-Landau band plasmons is accompanied by SdH type of oscillatory behavior 



in 1/B . These oscillations [2j, 122J are not with constant period in 1/B (which ex- 
hibits significant effect for small value of B and corresponding large value of 1/B) due to 
confinement potential acting in the x-direction and also show the depopulation, and cross 
over effects [14] on magnetoplasmons from density modulated two-dimensional electron gas 
(DM2DEG) |23j] to density modulated one-dimensional electron gas (DM1DEG). 

SdH type of oscillations result from the emptying out of electrons from successive Landau 
bands when they pass through the Fermi level as the magnetic field is increased. The 
amplitude of the SdH type of oscillations is a monotonic function of magnetic field, when 
the Landau bandwidth is independent of the band index n. In DM1DEG considered here, 
the Landau bandwidths oscillate as a function of the band index n. It is to be expected 
that such oscillating bandwidths would effect the plasmon spectrum of the intra-Landau 
band type, resulting in another type of oscillation. These oscillations are not with constant 
period in 1/B, because at small value of B and corresponding large value of 1/B cyclotron 
diameter exceeds the characteristic length of the confining potential. 

For the excitation spectrum, we need to numerically solve equation (15) for all vectors, 
frequencies, magnetic field, and confinement potential. We will consider the case of weak 
modulation (Vo/Ep << 1) and long wave length. In these limits we can solve equation (16) 
analytically for zero temperature. We expand the coefficient C nn i ( ) to lowest order in 
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its argument with the result 



i 2ne2 - 1 f mn °^ff \\ art 

u n 

The term in parentheses is easily recognized as the unmodulated particle density no = 
^^-'Yl l f{sn)i where n is sum over all occupied Landau bands. Defining the plasma frequency 

n 

through ui 2 D = 2 ™m ^' we fi nan y obtain the inter-Landau-band plasmon dispersion relation 

U! 2 

1 = ^2^2 or oj 2 = + oj 2 Di which is the ordinary one-dimensional plasmon dispersion 
relation. 

The intra-Landau-band plasmon dispersion relation for zero temperature reduces to 1 = 

~2 

where 

" = sin2 (^-o)x £ I ^ I V^ASf (1 - A n ), (18) 

with A n = - Af6>(l - A n ), A n =| |,0(x) the Heaviside unit step function. 

We have derived the expression for u (equation 18) under the condition oo »\ e(n, xq + 
x'q) — e(n,x ) | as x' — > which leads to a relation between the frequency and the Landau 
level broadening u »\ 2V n sin(^Xg) sin[(^)(x + ^)] |. This ensures that ImU (q,u) = 
and the intra-Landau-band magnetoplasmons are undamped. For a given V n , this can be 
achieved with a small but nonzero q y (recall that x' = — ) ■ 

In general, the inter- and intra-Landau-band modes are coupled for arbitrary magnetic 
field strengths. The general dispersion relation is : 



u 2 -tt 2 OJ 2 
This equation yields two modes which are given by : 

X (fig + C^D + - 2 _ 2 Q ^)}V2 

which reduces to: 

u 2 + = n 2 + 

and 

c<j 2 = <jj 

with corrections of order u /VLq and /w^ . So long as | y n |< Q , mixing of the inter-and 
intra-band modes is small. Only the intra-Landau-band mode ( a;) will be excited in the 



frequency regime Qq > & ~ I V n \ ■ We now present the results for the oscillatory bahavior of 
the intra-Landau band magnetoplasmons as a function of 1/B and the confinement energy. 

IV. NUMERICAL RESULTS AND DISCUSSION 

The intra-Landau-band plasma frequency given by equation (18) is shown graphically 
in Figure (ljas a function of 1/B for two different values of confinement energy, using 
parameters [3, S, Q]: m* = 0.07m e , k = 13.6, % = 6 x 10 15 m -2 , a = 500nm, and Vq 
= 1.0 meV; also we take q x = and q y = 0.01/cp, with hp = (27TO.D) 1 / 2 being the Fermi 
wave number of the unmodulated 1DEG in the absence of magnetic field. The modulation 
induced oscillations are apparent, superimposed on SdH-type oscillations. These oscillations 
are not with constant period in 1/B (which exhibits significant effect for small value of B 
and corresponding large value of 1/B) due to confinement potential acting in the x-direction. 
They have longer period and much reduced amplitude. It must be noted that this result is 
obtained in the regime of weak modulation and long wavelength. Another condition for the 
observation of these oscillations is that the coupling between inter-Landau band mode and 
intra-Landau band mode must be small. These modes involve different energy scales uj > Qq 
for the former and u ~| V n \< Qq for the latter. Mixing of these modes can be minimized 
by controlling the degree of density modulation and by applying an appropriate magnetic 
field and confinement potential. 

The origin of two types of oscillations can be understood by a closer analytic examination 
of equation (18). In the regime, f2 > \V n \, the unit step function vanishes for all but the 
highest occupied Landau band, corresponding , say , to the band index N. The sum over n 
is trivial, and plasma frequency is given as u =\ Vn \ 1 ^ 2 (1 — A 2 N ) 1 ^0{1 — An). The analytic 
structure primarily responsible for the SdH type of oscillations is the function 6(1 — An), 
which jumps periodically from zero (when the Fermi level is above the highest occupied 
Landau band) to unity (when the Fermi level is contained with in the highest occupied 
Landau band). On the other hand, the periodic modulation of the amplitude of the SdH 
type oscillations shown in Figures (1 & 2) is due to the oscillatory nature of the factor 
| Vn I 1 / 2 , which has been shown in a two-dimensional system to exhibit commensurability 



oscillations 
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241 ] . In our case, these oscillations are not with constant period in 



1/B. This clearly indicates the one dimensional character of our theory. In a DM2DEG 
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the number of occupied Landau levels increases with decreasing B, leading, ideally, to an 
infinite number of SdH type of oscillations periodic in 1/B 23j . In a 1DEG system however, 
only a finite number of ID subbands are occupied at B = 0, giving rise to finite number of 
SdH type of oscillations and deviations from the 1/B period, because with increasing B the 
D density of states increases and the hybrid ID subband Landau levels are depopulated 
4], 13, 2s| - In the extreme 2D regime (uq << uj c ), the Fermi energy goes to the bottom of 



the ID Landau subband. If we lower the confinement potential the magnetic confinement 
overcomes the potential confinement, hence we are in the original 2D regime 23j. On the 
other hand if we increase the confinement potential, the confinement potential overcomes 
the magnetic confinement and we have a crossover 14| from a two-dimensional system to 
a one-dimensional system. In figures (1, 2) we have plotted the intra-Landau band plasma 
frequency as a function of the inverse magnetic field for two different confinement energies 
given by Eq. (18). 



V. CONCLUSION 



We have determined the intra-Landau band plasmon frequency for a density modulated 
quasi-one dimensional electron gas in the presence of a magnetic field employing the SCF 
approach . Furthermore, we have seen the oscillations of the intra-Landau band plasma 
frequency in ID regime as a function of B~ l , their origin lies in the interplay of the three 
physical length scales of the system i.e. the modulation period, confinement length and 
cyclotron diameter at the Fermi level. When a strong magnetic field is applied, our model 
recovers complete Landau quantization and for a very high magnetic field our results are 
comparable with extreme 2D regime (f2 << uj c ). 
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